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This is the third of a series of papers dealing with factors affecting the accuracy of 
measurment of lens distortion by various visual and photographic processes. This paper 
presents the results of a study of the visual modified goniometric method which employs a 
viewing theodolite and accurate linear scale. Sources of error are investigated and discussed. 
Measurements are reported on the same lens that was used in the analysis of two visual 
and one photographic distortion measuring processes. A comparison of results by the four 
different methods is given. It is concluded that accurate values of lens distortion can be 
obtained by the modified goniometric method provided care is taken to use centered aper- 
tures and to preserve known angular relationships between telescope pointings and successive 
locations. 



1. Introduction 

The accurate evaluation of radial distortion in 
the focal plane of airplane camera lenses has been 
of i prime importance since the advent of mapping 
programs that make extensive use of photographs 
made with airplane cameras. Radial distortion, as 
herein used, is the linear displacement of an imaged 
point from the point where it would have been 
imaged if the lens were perfectly orthoscopic. This 
displacement is measured in the image plane of the 
lens along a radius drawn from the intersection point 
of lens axis and image plane. The presence of dis- 
tortion impairs the usefulness of the photographs in 
precision mapping. For this reason, specifications 
have been generally established which set limits on 
the magnitude of the permissible distortion for 
lenses used in various types of aerial photography. 
To determine compliance with these specifications, 
the values of the distortion are measured for each 
lens. Various methods have been developed for 
making these measurements and it is important 
that these produce reliable results. 

A study of several methods commonly used in 
the evaluation of lens distortion has been in progress 
at this Bureau for some time. The methods 
investigated thus far are: 

A. Photographic-precision lens testing camera. 

B. Visual — nodal slide optical bench with col- 

limator. 

D. Visual — inverse nodal slide bench with point 
or line target in rear focal plane of lens. 

F. Visual — modified goniometric with linear scale 
in rear focal plane of lens. 



In this study, measurements were made on the 
same lens by all four methods; the results were 
compared and attempts were made to locate the 
sources of error most likely to impair the reliability 
of results obtained by each method. The results 
of these investigations have been reported for 
methods A, B, and D; 1,2 the present paper gives 
the results for method F. 

A comparison of the values of the distortion 
obtained by methods A, B, D, and F are given in 
table 1 and figure 1. The values shown for A, B, 
and D are the final accepted values for these methods 
a It er all known sources of error have been eliminated 
or minimized; the values shown for F are the results 
of a single determination before any attempt was 
made to reduce or eliminate nonobvious sources 
of error. 

Examination of the values of distortion shown in 
part (a) of table 1 shows that while there is good 
agreement among corresponding values obtained by 
A, B, and D, the values obtained by method F 
depart markedly from those yielded by the other 
three. This is shown more clearly in part (b) of 
table 1 where the departures for each method from 
the average of values Z> 3 , obtained by A, B, and D 
are given. The values for F were not included in 
the average as it was believed that this set of values 
was affected by errors markedly in excess of the 
random errors of the order of ± 2\x that affected the 
others. This is further justified by considering the 
values of the departures from the average. Here the 
maximum departure of any value from the average 



* F. E. Washer, W. P. Tayman, and W. R. Darling, Evaluation of lens dis 
tortion by visual and photographic methods, J. Research NBS 61, 509 (1958) 
RP2920. 

2 F. E. Washer and W. R. Darling, Evaluation of lens distortion by the inverse- 
nodal slide, J. Research NBS 63C, 105 (1959). 
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Table 1. Measured values of the distortion in microns, referred 
to the calibrated focal length for wide angle lens No. 3 by four 
different methods. 

The values for A, B, and D in part (a) are final accepted values while those for 
F are for a single preliminary determination. The nominal focal length of the 
lens is 152 mm. Departures from the average, AD=A-Da, IH — Bb, etc., are 
shown in part (b). 
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Figuke 1. Variation of the departures, AD, from the average 
value of distortion, D3, with angular separation fi from the 
axis. 

The variation for each of the four methods are shown in the frames designated 
A, B, D and F. 



is ±5/z for A, B, and D, while departures ranging 
from —29 to + 17/z are present for method F. This 
is shown still more strikingly in figure 1, where it is 
clear that any random or systematic errors do not 
cause the values for A, B, and D to depart more than 
± 5/z from the average and such departures as exist 
are for the most part negligible ; it is also clear from 
figure 1, that the values obtained by method F are 
affected by large errors somewhat cyclic in nature. 
Because of the wide variance shown in the pre- 
liminary analysis, it appeared worthwhile to make a 
careful study of method F to determine the causes 
of these discrepancies. 

2. Study of the Modified Gonic metric 
Method (Method F) 

2.1. Theory 

In the usual goniometric method a device similar 
to a large spectrometer is used. The lens in its 
camera is so mounted that its front nodal point lies 
in the axis of rotation of the observing telescope. 
A scale in the focal plane of the camera is observed 
through the telescope and settings made on successive 
scale lines having known separations. For each such 
setting, the circular scale of the spectrometer is read 
in terms of angle. From the known separation of 
successive scale lines and the corresponding angular 
observations, it is possible to determine the focal 
length of the lens and the values of the distortion 
at selected points. 

The modified goniometric method, as here used, 
is similar to the method described by Merritt. 3 
The relations from which values of the focal length 
and distortion are obtained are the same for both the 
usual and modified goniometric methods. If a 
graduated scale is located in the focal plane of the 
lens, and the separations of successive scale divisions 
from the central division assumed to be located on 
the lens axis are designated as a u a 2 , . . . a n , and the 
angular separations of the corresponding divisions 
as viewed through the front of the lens be designated 
as ft, j8 2 • • • 0«, then 



f 1=ai cot 0!. 



(1) 



where j x is the equivalent focal length of the lens. 
Using this value of j and known values of 0! the 
distances to the other scale divisions are computed 
and the distortion, D, is given by the relation 



D n =a n —f l tan f3 n . 
2.2. Description of the Apparatus 



(2) 



The arrangement of the apparatus is shown in fig- 
ures 2 and 3. The lens is mounted in the nodal slide 
located on one leg of the T-beneh. The scale, a cali- 
brated meter bar, is placed in the focal plane of the 



3 E. L. Merritt, Methods of field camera calibration, Photoerammetric Eng. 
XVII, 610 (1951). 
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Figure 2. A rrangement of apparatus for method F, the modified 
goniometric method. 

The scale at the right is located in the focal plane of the lens mounted on the 
nodal slide in the center of the picture. The theodolite on the left faces the front 
of the lens under test. For the position shown, settings can be made over the 
range /3 = -fll to -11 degrees. 




Figure 3. Same as figure 1 except theodolite has been displaced 
along the bench for making measurements in an extra axial 
region. 

lens and placed as nearly normal as possible to the 
lens axis which is initially made parallel to the bench 
ways. The Wild theodolite is mounted on a saddle 
located on the cross arm of the T-bench. It is so 
placed that its objective is as near as feasible to the 
front of the lens under test. For this particular ar- 



rangement, settings can be made from +30 to —30 
mm on the scale. For wide angles the arrangement 
shown in figure 3 is used. The appearance of the 
scale of the meter bar as seen through the observing 
telescope is shown in figure 4 for three angles of 
viewing. 

2.3. Results of Measurement 

It is obvious that an entire diameter of the picture 
area of the lens cannot be investigated for any one 
position of the theodolite. In order to cover all 
points along a single diameter, the theodolite must 
be traversed along the bench in a series of steps. 
This gives rise to alternate possibilities. The theod- 
olite may be held in a single position with respect 
to the bench while a series of observations are made 
at selected points along the scale and then moved 
for successive series or it may be moved between each 
setting. The first process seemed easier in the early 
part of the work and was used. The results of meas- 
urement taken by this process are shown in the 
following seel ion in table 2. 

The center of rotation of the observing telescope 
is separated by a very appreciable amount from a 
vertical line passing through the front nodal point 



t 



\ 



AXIS 




Figure 4. Scale of the meter bar as seen through the viewing 
telescope. 
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Table 2. Measured values in the object space of the angular 
separation, £, from the axis corresponding to the 10 mm inter- 
vals of a standard scale located in the focal plane of wide angle 
lens No. 3 

Values are shown for corresponding positions along the same diameter to the 
right and left of axis. The position of the theodolite for each series of settings 
is indicated by the bench scale readings. The change in the value of p that 
occurs on shifting of the theodolite position is shown under A/S. 
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of the lens. Consequently for a single location of 
the theodolite on its bench it is possible to make set- 
tings on a comparatively few divisions of the scale 
located in the focal plane of the lens under test before 
vignetting cuts off the view. This is indicated in 
table 2, where the successive readings of the bench 
scale show that the theodolite occupied seven differ- 
ent positions to cover an angular range of /?=46° 
to the right and left of axis. When the theodolite 
is moved from one position to another along its 
bench, it is customary to make the first setting of 
the new position on the same line of the meter bar 
as the last setting of the former position. This is 
indicated in table 2 by the double entries for meter 
bar divisions 30, 90, and 130. The angular values 
obtained for these double entries do not coincide but 
show differences of varying amounts which are tabu- 
lated under the caption A/3. If it is assumed that 
the theodolite always is moved parallel to itself, these 
differences A/3 should equal zero. The alternative 
explanations are: (1) Curvature of the bench ways 
and (2) aberrations arising from an asymmetrical use 
of aperture. If the second cause is regarded as neg- 
ligible, then subtracting the decrement, A/3, from the 
observed values of /3 should eliminate any errors aris- 
ing from (1). This was done and the values of @ so 
determined are given in table 3. This table also 
shows the values of equivalent focal length and dis- 
tortion, D R , to right of axis and D L to the left of axis. 
The values of distortion in table 3 are not directly 
comparable in the form shown to those given in 
table 1. Consequently it was necessary to evaluate 
D R and D L for values of /3 at 5° intervals. This was 
done by linear interpolation and the values are 



Table 3. Values of the distortion, Dr and Dl, referred to the 
equivalent focal length (f= 152.518 mm), based on the meas- 
urements shown in table 2 

Values of /S are adjusted by adding the measured values of A/S to the measured 
values of /S for each bench setting. 
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Table 4. Values of the distortion, De and Dc for wide angle 
lens No. 3 for values of /3 varying in 5° steps 

These values are derived from table 3. 
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shown in table 4. The asymmetry of the distortion 
pattern results from a small tilting of the scalebar 
with respect to the focal plane of the lens. This 
effect can be readily compensated at the proper 
stage by computation. 4 The value of the distortion, 
D E) referred to the equivalent focal length is given 
by the relation 



D E 



D R +D L 



(3) 



Values of D' Cj the distortion referred to the cali- 
brated focal length are derived from D E in the usual 
manner. Finally these values are corrected with 
the aid of the equation 



D C = D' C -AD T , 



\4) 



where AD T is the correction for the tilt of the scale 
bar and D c is the value of the distortion referred to 



4 F. E. Washer, Sources of error in various methods of airplane earner (cali- 
bration, Photogrammetric Eng., XXII, 727 (1956). 
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the calibrated focal length which is accepted as the 
best value obtainable by this method of analysis. 
Comparison of the values of D c shown in table 4 
with the values of the average D 3 for methods A, B, 
and D given in table 1 shows marked disagreement. 
The disparity between corresponding values is so 
great that it is apparent that the measurements are 
being affected by some uncompensated source of 
error. 

It seems probable that the aberrations arising 
from asymmetrical use of apertures contribute 
heavily to the observed disparities. With this in 
mind it seemed worthwhile to neglect errors arising 
from bench curvature and to attempt to reduce the 
aberration effects by averaging the value of £ in the 
double-entry regions of table 3 on the hypothesis 
that the two separate values were obtained for 
opposite conditions of vignetting and that the aver- 
age might approximate the value obtained for cen- 
tered apertures. Tins was accordingly done and 
the values of /3 so determined are listed in table 5. 
No adjustment was made for those values of p at 
the non transition points. The values of distortion 
are computed and are shown in table 5. D R and D L 
are evaluated at 5° intervals and are tabulated in 
table 6, together with D E and D c . The values of 
D c obtained in this manner are nearer to the values 
shown in table 1 than are the values given in table 4 
which indicates that this method of reducing the 
data is superior to the method used in deriving 
table 4. The disparities that remain are due in 
part to some noncompensated bench curvature and 
in part to observation arising from unsymmetrical 
use of apertures. On making an adjustment, AD, 
for the error arising from the measured departure of 



Table 5. Values of the distortion, Dr and Dl, referred to the 
equivalent focal length (f= 152.518 mm) based on the meas- 
urements shown in table 2 

No adjustment if made for the values of 8 except to average the two values of 
6 at each transition point. 
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Table 6. Values of the distortion, Dr and Dl for wide 
angle lens No. 3 for values of varying in 5° steps 

The values of Dr (distortion referred to the equivalent focal length) and Dc 
(distortion referred to the calibrated focal length) are also given together with the 
adjustment of pointing error aDr and the final vnlue of distortion for the method, 
Df. These values are derived from table 5. 
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the bench from straightness, the value D F represents 
a value free from error arising from bench curvature. 
This value, D F , is still somewhat different from the 
average given in table 1 showing that the aperture 
error is appreciable. 

2.4. Errors Arising From Asymmetric Apertures 

In an earlier paper (see footnote 3), an extended 
discussion of the effect of asymmetrical apertures is 
given. It is evident from a consideration of the 
arrangement of apparatus used in the modified goni- 
ometer method that the apertures are being used 
asymmetrically for most settings. To show that 
this is so, a series of photographs of the effective 
aperture of the lens under test as seen in the exit 
pupil of the viewing telescope was made. Some of 
these photographs are reproduced in figure 5. The 
dark circles indicate the exit pupil of the telescope; 
the light circles or ovals show the effective aperture 
of the lens under test as seen in the exit pupil; and 
the irregular light areas are caused by random reflec- 
tions of stray light hitting the interior walls of the 
telescope tube. The principal interest in these 
photographs is to note the different relative positions 
of the lens aperture as seen in the exit pupil. Refer- 
ring to table 1, it may be noted that for bench scale 
reading 569 mm, and for meter bar reading of 30 mm, 
the aperture is at the right side of the pupil and 
partially obscured while for the second entry of 30 
mm at bench scale reading of 639 mm, the aperture 
is at the left of the pupil and cut even more. The 
difference in /3 for these settings is 0.0280°. Similar 
conditions exist for meter bar readings of 90 and 130. 
It is obvious that in proceeding from meter bar read- 
ing +30 to —30 the lens aperture is centered for 
only one setting at mm, and an asymmetric con- 
dition exists at all other positions in the range from 
+ 30 to —30. One would expect that values of the 
distortion obtained under these conditions would 
show large random variations from the true value. 
This is indicated in table 7 which shows a comparison 
between values of D 3 from table 1 and method F 
from table 6. 
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Figure 5. Photographs showing the relative position of the effec- 
tive aperture of the lens under test in the exit pupil of the 
measuring telescope for a series of bench and meter bar readings. 



Table 7. Comparison of the values of distortion Df obtained 
by method F from table 6 with the average values D3 from 
table 1 

For these values of T>f, the exit apertures of the lens and entrance aoerture of 
the" telescope are centered for values of /9 near 0°, 23°, 36°, and 44°. All values 
are in microns. 
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2.5. Method of Eliminating Effect of Unsymmetrical 
Apertures 

The improvement in reliability of results shown 
in table 6 over those shown in table 4 appeared to 
warrant additional measurements using a method 
calculated to minimize any effects arising from 
unsymmetrical use of apertures. Accordingly new 
measurements were made which are given in table 



8. In these measurements, it is assumed that the 
movement along the crossbench does not alter the 
reliability of the angle measurements, so the the- 
odolite was displaced for each setting in such a 
manner that the exit aperture of the lens under test 
appeared centered in the exit pupil of the telescope. 
The appearance of the exit pupil is shown in the 
right hand column for a few cases in figure 5. Each 
setting was repeated 5 times so as to determine the 
probable error of the mean for each value of 0. 
The probable error of the mean for the values of 
was found to be approximately ±2 seconds which 
corresponds to a maximum error in the distortion 
values of ±0.003 mm. 

The values of the distortion, D R and D Lj were 
determined from the measured values of ft cor- 
responding to the scale divisions of the meter bar 
and are given in table 8. These values were used 
to determine the values of D E (distortion referred to 
the equivalent focal length) at 5° intervals which 
are given in table 9. This table also shows the 
values of D c (distortion referred to the calibrated 
focal length). While the error arising from bench 
curvature was assumed to be negligible, it is known 
that this is not quite true. However, the magnitude 
of this error lias been measured for each position 
occupied by the theodolite and an appropriate cor- 
rection is made under the caption AD E (compensa- 
tion for pointing error). Finally the value of D F is 
shown in the last column. 

Comparison of the values of D'p from table 9 
with that of D z from table 1 is shown in table 10 
together with the difference Dz—D' F . Consider- 
ation of the values of D 3 —D' F leads to the belief 
that systematic error of some sort is still present. 
A possible source is curvature of the meter bar. 

Table 8. Measured values in the object space of the angular 
separation, /3, from the axis corresponding to 10 mm intervals 
of a standard scale located in the focal plane of wide angle 
lens No. 3 

Conditions the same as for table 2, except apertures are centered for each setting. 
Each value of is the average of 5 settings; the average probable error of ihe 
mean is ±0.0056° or 2.0 seconds which corresponds to a maximum distortion 
error of ±0.003 mm. The values of the distortion are also shown. 
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Table 9. Values of the distortion Dr, Dl, and De derived 
from table 8 

The values of Da, the distortion referred to the calibrated focal length and 
the values of D'f derived from Dc after compensation for pointing error AD v. arc 
also shown. All values are in microns. 
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Table 10. — Comparison of the values of Dp from table with 
D :j from table 1 

The magnitude of the correction Dc required to compensate assumed curvature 
of target scale (sagitta=22 m )and final accepted value Dp for method F are also 

shown. All values are in microns. 
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A computation was made, to determine the probable 
magnitude of the assumed curvature (see footnote 4). 
A value of 22ju for the sagitta of the curve seemed 
satisfactory and reasonable and is used to determine 
the values of AZ) C . The value of I) F shows the 
effect of making- this correction and is accepted as 
the best set of valines of distortion obtained by 
method F. 

2.6. Review of Results Obtained by Four 
Different Methods 

In the course of this investigation, extended 
series of measurement of the values of radial distor- 
tion have been made on a single lens using four 
distinctly different methods designated A, B, D, and 
F. Table 11 shows, in section (a), the final value 
accepted for each method together with the average 
for the four. Part (b) of table 11 shows the depar- 
ture from the average for each of the four methods 
together with the probable error of a single deter- 
mination for a given method with respect to the 
average for all methods. The departure from the 
average for each of the four methods is also shown 
in figure 6. In the course of the investigation, it 
has been shown that the PE S for values obtained 
by any one of the methods is approximately ±3 p. 
It is shown in table 11, that the departure of the 



Table 11. Pleasured values of the distortion in microns referred 
to the calibrated focal length for wide angle lens No. 3 for four 
different methods 

The values for A, B, D, and F in part A are the final accepted values obtained 
by each of the four methods. Departures from the average D for the four are 
shown in pari (b) together with the value of PE» t the probable error for values 
obtained by a single method. 
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(b) Departures from the average for method 
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Figure 6. Variation of the departures, AD. from the average 
value of distortion with angular separation j3 from the axis. 

The values of AD arc computed witli respect to the average value of D for the 
four methods A, B, D, and F using the final accepted values for method F. 

values obtained by a single method from the grand 
average for four methods does not usually exceed 
±3 p. 

It can therefore be concluded that the results 
obtained by any one of these methods are likely to 
be as good as those obtained by any other of the 
remaining three methods, and moreover the errors 
in the values of distortion are unlikely to depart 
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by an amount in excess of ±5 fx from tlie true value. 
This statement is made on the assumption that 
proper care is taken in measurements for any of the 
methods and that all conditions are equivalent. 

3. Conclusions 

While the equivalence of the four methods studied 
can be regarded as established as a result of these 
investigations, there are other faetors that are worthy 
of mention in connection with them. It can be safely 
stated that method B is the easiest of the methods 
to use and requires the simplest apparatus, namely 
nodal slide, visual optical bench, and collimator. It 
is also relatively easy to understand. Method D is 
quite easy and simple but more difficult to compre- 
hend because of the inverse manner of the use of 
the nodal slide. In addition, unusual care is required 
to insure proper centering of apertures unless a large 



diameter viewing telescope is used. Method F while 
reliable is the most difficult of all. There appear to 
be more sources of error likely to creep in than for 
method B and D, and in addition, it requires more 
extensive use of trigonometric tables and interpola- 
tion than an}' of the other three methods. Method 
A is possibly the easiest to understand. It does how- 
ever, require complex equipment. It is slower than 
method B when there is only one lens to measure. 
However, as the number of lenses to be tested in- 
creases, it becomes more efficient as no time need be 
lost waiting for negatives to dry before measurement 
and it is possible to complete measurements on 10 
lenses quicker by method A than by method B. 
Method A has, of course, the major advantage that 
it approximates more nearly conditions of use than 
saiy of the other methods. 
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